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ABSTRACT 

Properties of stars undergoing pulsation such as the well known root-mean- 
density scaling relation can be useful when trying to match the observed 
properties of a particular star. It is often assumed that this relation is valid for 
p mode frequencies in rotating stars. To examine the change in frequency with 
rotation and mass, we have studied oscillation frequencies of two-dimensional 
uniformly rotating zero-age main sequence stellar models in the delta Scuti 
mass range. We identified axisymmetric p and g modes for non-rotating models 
and then traced them as the rotational velocity was increased. We considered 
a rotation sequence of ten models for four different masses, with the largest 
rotation rate being about 200 km s“^. The models were required to have the 
same surface shape between all masses for a given rotation rate. We find that 
scaling relationships exist among the oscillation frequencies of the same mode 
for different masses when the models have the same shape. For p modes, this 
scaling closely follows the period root-mean-density relation found in spherical 
stars. The g modes also scale between models of the same shape, with the 
scaling reflecting the change in properties outside the convective core as the 
stellar mass increases. These scaling relationships can be particularly useful in 
finding specific stellar models to match the oscillation frequencies of individual 
stars. We also find that the large separation scales approximately with the root 
mean density as the rotation rate increases, although the individual mode 
frequencies do not. 
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1 INTRODUCTION 


It has long been recognized that stellar oscillations provide a unique opportunity to study 


the interior structure of stars. Thanks to the space based missions MOST (Walker et ah 


2003), CoRoT (Baglin et ah 2001), and Kepler (Basri et al. 2005), we now have accurate 


frequencies for a large number of stars, including intermediate mass main sequence stars 


such as the 6 Scuti stars (e.g.. 

Michel et al. 2008; 

Kjeldsen et al. 

2010 

). While much effort 

has been expended modeling these stars (e.g.. 

Templeton et al. 

1997; 

Breger 

2000 

Suarez 


et al. 2014 Ouazzani et al. 2015), a global understanding of the fine structure of these stars 


still eludes us. 

A signihcant complication for many of these stars is rotation. Very small rotation re¬ 
moves the degeneracy of the non-axisymmetric modes, but only modest rotation is required 
before the mode splitting is comparable to the large separation, complicating any analy¬ 


sis of an observed frequency spectrum (Suarez et al. 2010 Deupree & Beslin 2010). Not 
much more rotation is required before the splitting of the non-axisymmetric modes becomes 
nonuniform. Rotation also alters the frequencies of the axisymmetric modes (e.g., Saio||1981 


Lignieres et al.||2006 Lovekin fc Deupree||2008). Finally, significant rotation deforms the star 


sufficiently so that the spectral energy distribution (SED) depends on the (usually unknown) 
inclination between the rotation axis and the observer, so that even placing the star in the 
HR diagram becomes nontrivial. The picture is brightened somewhat by the advent of op¬ 
tical interferometry, which allows us to determine both the inclination and, to some extent, 
the amount of rotation through the observed distortion of the stellar surface for a limited 

Belle et al.||2001 Domiciano de Souza et al.||2003 Aufdenberg et al. 


number of stars (e.g. 


van 


2006 

Monnier et al. 

2007 

Zhao et al.|2009 


For a rotating star the horizontal and radial variables of the perturbations in a linear 
pulsation analysis are no longer separable as they are for a non-rotating star. In effect this 
means that the latitudinal component of the eigenfunction can no longer be given by a single 
Legendre polynomial, but must be regarded as a sum of Legendre polynomials. This has led 


to two dimensional solutions to the eigenvalue problem (e.g., Clement 1998 Yoshida & 


Eriguchi 2001 Espinosa et al. 2004 Lignieres et al.||2006 ; Lovekin fc Deupree||2008 ; Lovekin 


et al. 2009; Ouazzani et al. 2012, 2015). 


Under the best of circumstances the task, then, is to find a rotating model which has 
the observed oblateness and oscillation frequencies and which, when viewed at the observed 












































































































































Scaling of oscillation frequencies in rotating stars 3 
inclination, produces the observed SED. Because the oblateness is tightly coupled to the 
mass, stage of evolution and the rotation rate, it is unreasonable to expect that just evolving 
rotating models from the ZAMS to the approximately right place in the HR diagram will 
satisfy all these conditions. Fortunately, we can simplify the process in various ways. 


Deupree (2011a) computed a number of ZAMS models with different mass and rotation 


rates and found that any two models which have the same surface shape (close to having 
the same value of H/Hcrit but not quite) have the same relative distribution of the effective 
temperature as a function of latitude. The effective temperature and effective gravity as 
a function of latitude and the surface shape are used to compute the SED that would be 


observed at a given inclination between the rotation axis and the observer (e.g., Slettebak 


et al.|1980 

Linnell & Hubeny|1994 

Fremat et al.|2005 

Aufdenberg et al.|2006i Lovekin et ah 


2006). This suggests that the deduced effective temperature and deduced luminosity, which 


one would obtain from the SED under the assumption that it originates from a spherical star. 


might have the same relative distribution as functions of the inclination as well. Castaneda 


& Deupree (2014) performed calculations of the SED and showed that this is the case as 


long as the models are not too dissimilar in structure (much like the period - root mean 
density relation applies). 

Such scaling allows the determination of the surface conditions of a rotating star with 
relatively little computation. The scaling algorithm assumes that the oblateness and incli¬ 
nation of the star is known from interferometric observations and that there is access to a 
rotating model with the observed oblateness. Besides the surface properties of this model, 
the (computed) deduced effective temperatures and luminosities as a function of inclination 
are required. Also needed are the deduced effective temperature and deduced luminosity at 
the observed inclination of the star, which can be computed from the observed SED in the 
same way as they were computed for the model. With this information, the ratio between 
the surface properties and the deduced properties of the model and the observed star can 
be calculated, leading to all the surface and intrinsic properties for the rotating star. This 
is equivalent to placing a non-rotating star in the HR diagram. See [Castaneda &: Deupree 


(2014) for more details regarding the scaling process. 

The next step is to compute a model that has the correct shape, surface effective temper¬ 
ature distribution, and total luminosity. Evolution tracks should get one close, and then one 
converges models with small perturbations to the rotation rate, total mass, and composition 
distribution to obtain the model. Of course, it is unlikely that this model will have the ob- 
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served oscillation frequencies, so we must now make corrections to the model to improve the 
agreement between the observed and model oscillation frequencies. This can be most easily 
performed if there are scaling relations between the oscillation frequencies of two models, 
much like the conditions for which the period-root mean density relation for non-rotating 
stars, to guide us. 


Previous studies have been able to relate properties of the oscillation spectra like the 


large separation with physical properties of the star, particularly the mean density (Reese 


et ah [2008 Garcia Hernandez et al.||2013 Garcia Hernandez et al.||2015) for stars close to the 


5-Scuti range of masses. This is can be an important tool when trying to hnd patterns that 


can aid the mode identihcation process. In particular, Lovekin et al. (2009) tested physically 
meaningful radii that would make the period-mean density relation, Q = P^Jm/R^, constant 
as the rotation rate increases. Using the radius at 40° seemed to keep Q constant as a 
function of the rotation rate, but its selection had no physically meaningful justihcation. 
For the models we are considering however, a natural consequence of comparing two stars 
with different masses but the same surface shape is that by dehnition their physical volume 
will be proportional between each other. Additionally, if the range of masses is small enough 
so that their internal structures are similar, this volume proportionality could potentially 
relate their oscillation properties directly, just like the period - root mean density scaling 
relationships for non-rotating stars. One purpose of this paper is to test the validity of this 
hypothesis. The mean density can be changed by at least three quantities associated with 
stars; the stellar mass, rotation, and evolution. With our models we can judge the relative 
importance of each in terms of scaling oscillation frequencies using the mean density 


Section 2 presents the procedure followed to obtain the oscillation frequencies of the 
models considered. In Section 3 we present an overview of the modes found in the frequency 
range considered, mainly highlighting some of the challenges encountered while each mode 
was followed as rotation rate increased. Section 4 presents the results of what happens when 
the oscillation frequencies between two rotating models that have the same surface shape are 
compared. Section 5 discusses the relative importance of the stellar quantities that can affect 
the mean density. Finally section 6 contains a review of the results, potential applications 
as well as the limitations of the relationships found. 
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Table 1. Rotating models considered 


Mass (M 0 ) 
1.875 2.000 2.250 

Shape (RpoieZ-^eq) Velocities (kms~^ 

2.500 

) 

^eff,pole/^eff,eq 

1.000 

0.0 

0.0 

0.0 

0.0 

1.000 

0.997 

35.0 

36.0 

36.0 

37.5 

1.003 

0.991 

62.0 

63.0 

65.0 

67.0 

1.009 

0.985 

83.0 

84.0 

87.0 

89.0 

1.015 

0.976 

105.0 

106.0 

109.0 

111.0 

1.024 

0.966 

125.0 

127.0 

131.0 

134.0 

1.034 

0.954 

146.0 

148.0 

152.0 

156.0 

1.047 

0.940 

165.0 

168.0 

173.0 

178.0 

1.064 

0.925 

187.0 

190.0 

195.0 

200.0 

1.079 

0.907 

207.0 

211.0 

217.0 

222.0 

1.098 


2 CALCULATION OF OSCILLATION FREQUENCIES 

We follow the same procedure for calculating the oscillation frequencies of rotating stars 


used by [Lovekin fc Deupreq ( |2008[ ) , [Lovekin et al.| ( |2009[ ) and [Deupree fc Beslin] ( |2010[ ) . 
The method is a two step process: First we use 2D rotating models calculated using the 


ROTORC code developed by Deupree (1990, 1995). The code solves conservation equations 
for mass, three components of momentum, energy and hydrogen abundance, as well as 
Poisson’s equation for the gravitational potential. The particular models included in this 


work are all a subset of the models presented in Deupree (2011a); they are are zero-age main 
sequence (ZAMS) models with uniform rotation and solar metallicity with masses between 
1.875 and 2.5 Mq, approximately covering the 5 Scuti mass range. For each mass we picked 
ten surface equatorial rotational velocities, keeping Rpoie/Req of the model constant between 
the different masses at every step in rotational velocity. A summary of the properties of these 
models is given in Table [T] It is notable that the ratios between the polar and equatorial 
effective temperatures are effectively the same for all masses for a given ratio of polar to 


equatorial surface radius (Deupree 2011a). The effective temperatures at each latitude are 
not determined by von Zeipel’s law, but rather are determined by a relation between the local 
effective temperature and local surface temperature. This makes d\ogg/dXogTew smaller, 
closer to 0.2 than to the value of 0.25 given by von Zeipel’s law. This smaller value is similar 


to that found in previous studies by various means (e.g., Monnier et ah 2007 Che et ah 
2011 Espinosa Lara & Rieutord 2011 Claret 2012[ ). 

We calculated pulsation frequencies using the linear, adiabatic, non-radial pulsation code 


developed by Clement (1998). This code expresses the eigenfunctions of each mode as a sum 
of up to eight spherical harmonics with a radial resolution of 500 zones. In particular, for 
this study we calculated modes with £ values up to 15 (even and odd parity modes can be 
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computed separately). The primary motivation for the number of basis functions is to obtain 
the frequencies of the low i modes accurately, not to compute the highest £ modes. For this 
reason, we choose to include in our analysis only values of i between 0 and 6. Similarly, 
we focus our discussion on the radial orders of p modes up to n = 8 in the non-rotating 
case because of our limited resolution near the model surface. We then follow these specihc 
modes as a function of rotation, although we do take into account avoided crossings between 
the modes we are interested in and the higher i modes. As was found by previous studies 


using this code (Lovekin & Deupree 2008 Deupree & Beslin 2010), we find that 8 spherical 
harmonics can yield accurate eigenfrequencies for selected modes and rotational velocities 
considered. Specihcally we found that the difference in computed frequencies of modes using 
6 and 8 basis functions is on average on the order of a few thousands of a percent for 
the lowest rotating models and about a tenth of a percent for the highest rotating models. 
These frequency differences indicate that at least for the modes selected, the number of basis 
functions is adequate. 

For rapidly rotating stars the process of classification of modes is more problematic than 
for non-rotating stars. The major reason is because the usual latitudinal quantum number, 
i, is no longer sufficient to specify the latitudinal variation of the eigenfunction. Similar to 


other studies, (e.g., Lignieres et ah 2006 Lovekin & Deupree 2008 Lignieres & Georgeot 


2009 Deupree & Beslin 2010) we address this problem by characterizing each mode with 


the value of i associated with the mode in the non-rotating model, £o. This tracing process, 
although useful, becomes progressively more difficult for more rapidly rotating models. For 
each model we identihed the same combination of g and p modes as well as several radial 
modes. In total, we computed 180 modes with even parity and 189 modes with odd parity. 
Table presents the non-rotating modes identified. Of these, the number of modes which 
match the analysis criteria described above is 78. Under the classihcation scheme developed 
by [Lignieres et ah (2006); Lignieres & Georgeot (2009); Reese et ah (2009) in which modes 
for rapidly rotating stars can be separated as either island, whispering gallery or chaotic 
modes, we expect our selected low degree modes to be mainly island modes. 


3 VARIATION OF MODE FREQUENCIES WITH ROTATION 

Significant research has been performed to determine the behavior of oscillation frequen¬ 
cies with rotation. Some studies have used the “traditional” approximation in which the 
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Table 2. Oscillation modes calculated for each model. The value of n denotes the number of radial nodes found in mode for 
the non-rotating case 


£ 

n-range {g mode) 

n-range (/ and p mode) 

0 


0-15 

1 

1 

1-16 

2 

1-2 

0-14 

3 

1-3 

0-15 

4 

1-4 

0-15 

5 

1-5 

0-15 

6 

1-7 

0-15 

7 

1-8 

0-14 

8 

1-9 

0-13 

9 

1-10 

0-13 

10 

1-11 

0-13 

11 

1-13 

0-13 

12 

1-13 

0-12 

13 

1-15 

0-12 

14 

1-16 

0-12 

15 

1-18 

0-12 


horizontal component of the angnlar momentum is neglected in the linearized momentum 


equations (e.g., Lee & Saio 1986, 1997 Townsend 2003). Two other prominent methods 
are perturbation techniques (e.g., Saio 1981 Dziembowski & Goode 1992 Souh et ah [1998 


Suarez et ah 2010) and, in principle the more reliable, numerical integrations with the hori¬ 


zontal variations assumed to be given by a sum of spherical harmonics (e.g., Clement 1998 


Lignieres et ah 

... 

Lovekin & Deupree 

2008 

|Ballot et al. 

2010) 


explored the problem using ray theory principles (e.g, Dintrans & Rieutord 2000; Lignieres 


& Georgeot 2009 Pasek et ah 2012 Prat et ah 2016) to compute oscillation frequencies in 


the asymptotic regime. Here we use Clement’s radial hnite difference integration approach 
in the radial direction. Furthermore, we will consider only axisymmetric modes. We shall 
distinguish different latitudinal modes by even though the modes are actually linear 
combinations of different £’s. When we use a specihc value of i, we intend this to mean the 
value to which a given mode could be traced back to at zero rotation. Figure shows the 
general evolution of oscillation frequencies with rotation. The p mode frequencies decrease, 
mostly because the stellar volume increases, and hence the mean density decreases, with 
rotation. The g modes, on the other hand, generally increase in frequency as the rotation 


rate increases (e.g.. Ballot et al.|2010), but by a much smaller amount. This hgure also shows 
a subset of low radial order g modes that increase in frequency with increasing rotation at 
low rotation, but then begin to decrease with frequency as the rotation increases further. 

To discuss these and other features, we divide the frequency range covered into four inter¬ 
vals: 1) low frequency which is dehned as having only g modes, 2) an intermediate frequency 
region which includes both g and p modes, 3) a higher frequency region which includes 









































































8 D. Castaneda and R. G. Deupree 

only p modes, but at frequencies sufficiently low that there are still sizable variations in the 
frequency difference between radial quantum numbers n and n+1 for different latitudinal 
modes, and 4) the high frequency regime in which the large separation is essentially the 
same for all latitudinal modes. We shall address each of these frequency regions in turn. 


The low frequency g modes in our set are perturbed much less by rotation than the p 
modes for the uniformly rotating models. This is because the g mode frequencies are most 
sensitive to the region just exterior to the convective core boundary while the rotational 
effects on the model structure are substantial only near the surface. Both the frequencies 
and the eigenfunctions in this interval are modestly perturbed, making mode identihcation 
relatively easy. The change in frequency from the non-rotating value for a given rotation 
rate is larger at lower frequencies, increasing by about 40% as the frequency is decreased by 
about a third. There seems to be little dependence with i. This means that avoided crossings 
are rare, but they are not entirely absent. It would appear that a necessary condition is that 
the g mode frequencies are already quite close in the non-rotating model. Interestingly, 


previous studies (e.g.. Ballot et al. 2012 Takata & Saio 2013) have found that some of 
these g modes with frequencies close to each other show particular rosette pattern in their 
internal structure. It is still not clear whether these modes are excited in real stars but we 
can conhrm that some rosette modes are part of our data set. Figure shows an example 
of one case found in our sample, in which a g mode with £ = 3 and n = 3. The frequency 
evolution of these modes, however, does not seem to differ from the non-rosette modes, 
making them only interesting in their internal structure. Finally, we note that even our 
lowest frequency modes are somewhat above the sub-inertial range (thereby limiting the 
impact of the Coriolis force) and that the difference between their pulsation frequencies and 
the corresponding non-rotating values increases as the square of the rotation rate, reflective 
of the centrifugal force. 


The region in which there are both p and g modes is the most complex. This region is 
relatively small in frequency width (the width increases as the rotation rate increases), but 
the mode density is relatively high and the different signs of frequency change for the p and 
g modes with increasing rotation make near resonances and avoided crossings comparatively 
common. We present an example of avoided crossings in Figure and also note from Figure 
[T] that the hrst p mode shows a change in its frequency evolution trend caused by one of 
these avoided crossings. Avoided crossings are relatively more common at higher £ because 









Scaling of oscillation frequencies in rotating stars 9 
higher ^ g modes for a given n have higher frequencies. Similarly, they are more common at 
lower n because lower n g modes have higher frequencies. 

It is also in this interval that we see the phenomenon mentioned above where a given 
g mode frequency increases with increasing rotation up to a point and then decreases with 
further rotation (see Fig. [^. These are mixed modes that have nodes in the interior where 
the g mode nodes are expected and also much nearer the surface where p mode nodes occur 
(e.g., Scuflaire|1974 Osakijl975 ). The frequency may actually begin to increase with increas¬ 
ing rotation before the node near the surface appears. This is similarly indicated in phase 
diagrams, such as dehned by Scuflaire (1974) and Osaki ( |1975 ), in which the latitudinal 
displacement or, equivalently, ^ + (the primes denote Eulerian perturbations) is plotted 
against the radial displacement. The latitudinal displacement and the pressure-gravity per¬ 
turbation are functionally equivalent in a non-rotating model, but not in a rotating one. We 
use the Eulerian pressure-gravity perturbation in part for convenience and in part because 
of the similarity of the latitudinal variation of this perturbation and the radial displacement. 
The phase diagrams are computed along radial lines (i.e., at a given latitude). At low rota¬ 
tion rates the phase diagrams at different latitudes show essentially the same clockwise or 
counterclockwise behavior, but the values of neither variable are the same. This is most no¬ 
ticeable at the surface, where the hnal point on the diagram may end in different quadrants 
for different latitudes. At higher rotation rates there may be different numbers of nodes at 
different latitudes, and in the phase diagram this may appear as loops which do not quite 
cross one or the other axes at some latitudes but do at other latitudes or as a direction 
reversal which crosses an axis for some latitudes but stops short of crossing at others. The 
relationship between the number of nodes and a specific mode identihcation becomes more 
murky at the highest rotation rates we use, although the motion in the phase diagram is 
still clockwise for g mode behavior and counterclockwise for p mode behavior with the type 
of mode behavior determined by where the nodes are in radius. The change of the phase 
diagram with rotation for a mode that will become mixed as the rotation rate increases is as 
follows. As the radius increases, the curve proceeds in a clockwise direction for g modes, but 
for these g modes that show a frequency decrease with increasing rotation, the motion be¬ 
comes counterclockwise at sufficiently large fractional radius, although it does not continue 
far enough to produce a node near the surface. Because the effect increases as the rotation 
rate increases, one expects that p nodes will appear at yet more rapid rotation. We do see f 
modes that are mixed in these non-rotating ZAMS models, and adding rotation eventually 
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leads to eigenfunctions in which the node pattern is more pronounced. Rotation also extends 
mixed modes to some low radial order g modes. There may be an avoided crossing before 
this mixed original g mode occurs. 

Avoided crossings also occur in the third frequency region where there are no g modes. 
The usual case is a low latitudinal order p mode decreasing its frequency with increasing 
rotation reaching the frequency of a higher latitudinal order p mode that is decreasing 
its frequency more slowly with increasing rotation. Some easily identihable crossing in our 
models appear to involve £ = 0, 4, 8, although other cases certainly occur, in agreement 


with the crossing pattern found by Lignieres et ah (2006). We show one such progression in 
Figure The internal structure illustrations on the right show how it is still fairly easy to 
differentiate between the two modes involved after the interaction. 

The hnal frequency interval is where n is sufficiently high that the large separation 
is virtually the same for all £ values. Thus, one has as many modes as one has spherical 
harmonics within a frequency interval given by the large separation. The large separation 


decreases very slightly with increasing rotation, as was found by Lignieres et ah (2006) and 


by Lovekin & Deupree (2008). Another feature is that the eigenfunction for a given I no 
longer changes very much from one value of n to the next for a given latitudinal mode. It was 
hoped that this frequency regime would be useful in mode identification, but the placement 
of modes within the large separation, including the ordering, may be different from one 
rotation rate to the next. For higher frequency modes techniques like the application of ray 
tracing methods have been successful at finding insights into the frequency patterns that 


Lignieres & Georgeot 

2009; 

Pasek et al.f 

2012) 


methods, however, do not replace the need for highly resolved two dimensional calculations 
that provide more exact eigenfrequencies and eigenfunctions. 

It is interesting to note that the trends for these intermediate mass ZAMS models are 


similar to those of the 10 Mq model of Lovekin & Deupree (2008) and the rotating polytropic 


models of Lignieres et ah (2006). 


4 VARIATION OF MODE FREQUENCIES WITH MASS 

The period-root mean density relation has long been a tool for scaling the oscillation fre¬ 
quencies or, alternatively, estimating the mean density for a model or star for which there is 
an object sufficiently similar in properties to make the comparison. It has often been tacitly 
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Figure 1. Rotational evolution of = 4 mode for M = 2.5Mq. Frequencies are scaled by the factor {AnGp^ei)^^'^<> where 
= Igcrn”^, to keep their values at around unity across models and Vj-ef = 100 km s“^. Some of the changes in direction, 
most evidently for the first p mode, are caused by interactions with other i modes as well as changes of properties of each mode 
as rotation is higher. Radial order increases to the left for the p modes and to the right for the g modes 


assumed to remain useful for rotating stars. Strictly speaking, one can only expect it to be 
applicable to p modes, and our discussion about the different behavior of p mode and g mode 
frequencies with rotation suggests that the period-root mean density law does not apply for 
g modes. Here we will determine the applicability of the period-root mean density law for p 
modes of rotating stars and examine the scaling of g modes. 

Figure [^presents the frequency ratio of the same modes for the ZAMS models of 2 and 
2.25 Mq. Figure]^ (left) shows the results for non-rotating models and Figure]^ (right) shows 
the same calculations for moderately rapidly rotating models. In each case the ratio of the 
two root mean densities is shown by the horizontal line. It should be noted that the two 
rotating models having the same surface shape means that the scaling of the stellar volumes 
will be exact and given by the cube of the ratio of the equatorial radii. It is clear that the 
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Figure 2. Meridional cross section of the pressure perturbation of a ^ mode (^ = 3, n = —3) showing the rosette-like structure. 
A colour version of this figure is available online 


effect of rotation on the p mode ratios is rather small and that the period-root mean density 
is still a useful indicator of how to change the model to obtain desired period changes. 

It is also clear that the g modes do not follow this rule, although the ratios also do not 
depend very much on rotation. We can perhaps obtain some insight into how this g mode 


ratio comes about by applying equation equation (A13) from Tassoul (1980) to the stable 
g modes with radial nodes located between the outer boundary of the convective core and 
the surface ionization regions 

^ w+yryjyyi 

^ {n + l/ 2 )'K - 6 i- 62 

Here N is the Brunt-Vaisala frequency, and the thetas are phase shifts which occur if either 
N or the density or their derivatives are discontinuous at the limits of the integration. 
Because N is taken to be identically zero in the convective core, its hrst derivative may 
be discontinuous at the convective core boundary. The integral is taken between two radial 
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Figure 3. Rotational evolution of £’s with even parity calculated for the 1.875Mq models. Vref Pref defined in Figure 

1. We have labeled each line with its corresponding “initial” i. There are a few cases of interaction between p modes whose 
frequency decreases as rotation increases and g modes whose frequency evolution behavior is the opposite. It is important to 
note that the lines do not cross, they correspond to avoided crossings. Modes in this region are the hardest to follow because in 
some cases modes can interact with each other two or three times before reaching the most rapidly rotating models included. 


locations which are consecutive zeros of iV, in this case the convective core boundary and the 
lower boundary of the second helium ionization region. Assuming that we have successfully 
identihed the same mode and the phase shifts are the same for each of the two models, the 
ratio of the two mode frequencies will be 

rr2a j\j dr 
Jrla r 




, ,, rrzo ]\Tdr 

Jrlb A- 

Performing the integrations with the two non-rotating models in Figure we obtain a 
ratio of 1.088, clearly close to the upper limit for the ratios of the linear, adiabatic frequencies. 
For all masses, the largest contributions to the integrals come from the region exterior, but 
relatively close, to the convective core. This is true both because N is relatively large and r 
comparatively small there. However, the contributions all the way out to the bottom of the 
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Figure 4. Mode interaction between modes with = 0 and £q = A for models of mass M = 2.25M0. Vj-ef pref 
defined in Figure 1. It is possible to see the change in “label” required to match the exchange of surface latitudinal variation 
for each mode. The 2D contour plots on the right show the before and after of the meridional cross section of the pressure 
perturbations of the interacting modes enclosed in the box. In the case of the pressure perturbations, the lower frequency mode 
is on the left. 


second helium ionization region are not entirely negligible, even if less prominent than the 
deeper regions. 

We also found that as rotation rate increased some of most interesting individual mode 
ratios were found in the frequency region that contains both p and g modes. For the non¬ 
rotating models, the modes with erratic ratio behavior with rotation are either the / mode 
or low radial order g modes, with no particular preference in latitudinal order. However, 
mode ratios in this region become signihcantly more complex with rotation. The ratios of 
some of these modes change very little with rotation, while the ratios for others change 
signihcantly. While no one feature accounts for the modes whose ratio changes signihcantly, 
modes bumped, modes with near resonances, and mixed modes appear to be more likely 
to display this behavior than others. The order at which modes appear in the non-rotating 
case in this region also seems to be diherent for every mass and can can also account for 
the erratic frequency ratio behavior. For the non-rotating cases about 12-18 modes out of 
the about 140 modes in this region, including all i values, change places. Most of the modes 
involved in this ordering change, however, are modes with i > 8 and not taken into account 
for the ratio analysis. 

We also see from Figure that the ratios do not change signihcantly as rotation is 
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Figure 5. Ratio between frequencies of stellar models with the same shape. Left presents the non-rotating case. Right shows 
the same ratio for two rapidly rotating models. The black horizontal line indicates the ratio between the root-mean densities 
of both models. The cross shown on the figure of the right indicates what the ratio between modes that have undergone an 
avoided crossing for one mass but not the other. See the text for more details. 


included (except as we have indicated), at least up to the rotation rates presented here. We 
suggest that one may use the period-root mean density relationship for p modes in rotating 
models with the same reliability and constraints as one uses it in non-rotating models. The 
lower frequency g modes also scale predictably; it is only the modes in this intermediate 
region which are less predictable. 

Finally note that in Figure we have followed the mode frequency evolution, not the 
shape of the latitudinal variation of the eigenfunction. This makes no difference in terms 
of the computed ratio unless one has the case where a particular mode has undergone an 
avoided crossing for one mass but not the other. These are rare occurrences. Following the 
shape in this case changes the frequency ratio for both modes. We show the approximate 
shift by the cross in Figure]^ (right). In terms of increasing rotation rate, the ratios look 
normal until we reach the rotation rate where the avoided crossing has occurred in one model 
but not the other. Once the avoided crossing has occurred in both models the ratios look 
normal again. 

Figure|^shows an overview of how the scaling of frequencies look for all masses considered 
using the eighth most rapid rotating models in our set with velocities around 170 km s“^. 
It can be seen how the g modes scale differently than the p modes in every case. To give a 
general idea of how the ratios of the p modes scale for the different mass combinations as 
well as to see how each ^ behaves, the average ratio between each acoustic mode with n > 1 
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Figure 6. General scaling behavior of the frequency ratios between the computed frequencies of the M=1.875 Mq and the 
other masses considered for the eighth most rapid rotating models in our set with velocities around 170 km s“^. Modes with 
£ values from 0 through 6 were included. It can be seen how g and p modes scale differently. The horizontal lines correspond 
to the theoretical ratio of the root-mean density relation for each pair of masses. The p modes fall close to these lines in every 
case, but depart from the lines progressively more as the relative difference between the two masses increases. 


and the theoretical root-mean density ratio was calculated. Most ratios fall within 0.1% and 
1% from the theoretical ratio, where the larger differences correspond to the ratios between 
models with the biggest mass difference. 


5 DEPENDENCIES OF THE Az/ - 


Several studies have been undertaken to examine the utility of the Au — s/p relation for 6 


Scuti variables (e.g., Reese et ah 2008; Suarez et ah 2014 Garcia Hernandez et ah 2015). 
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Here Av is taken to refer to frequency differences with An = 1 and A^ = 0 for relatively 
small n instead of in the limit of large n. We discuss the relative importance of three factors 
that affect this relationship: evolutionary state, rotation, and mass for models with the same 
rotational shape. 

The largest variation in the mean density comes from the evolution of the star along 
the main sequence, dropping by a little less than a factor of six from the ZAMS to the 
end of core hydrogen burning for a 2Mq non-rotating model. In contrast, the change in 
mean density between our 2Mq non-rotating model and a 2Mq model rotating near critical 
rotation ( Deupre^|2011a ) is less than a factor of two. The change between the mean density 
of the non-rotating model and the highest rotation rate we consider in this paper is only 
about 10%. We obtain the same mean density as that for our most rapidly rotating model for 
an evolved model with Xc = 0.64 (the ZAMS models have Xc = 0.70). Finally, we note that 
the mean density changes by about 17% between our 1.875Mq and 2.5Mq ZAMS models 
(the ratio is the same for all rotation rates because of the way the models were constructed). 

We compare the relation between Av and y/p as a function of mass for the same evolution 


state, rotation, and the relation found by Garcia Hernandez et al. (2015) and by Suarez et al. 


(2014) in Figure]^ Like Garcia Hernandez et ah, we computed Av at sufficiently high n so 
that we could avoid our second frequency region where the p and g modes overlap and the 
density of modes is comparatively high. This was done by locating the highest frequency 
mode with g mode characteristics in our most rapidly rotating models and considering the 
next six modes for each i that can be traced back to modes with 0 ^ ^ 5 in the non-rotating 

model. This gives us the modes with 3 ^ n ^ 8 for which we compute Az/ averaged over 
the hve large separations for each of the six values of i. The modes we use have sufficiently 
large n that the variation of Az/ with n and ^ is relatively small in the narrow range we use. 
We do note that, for lower radial order p modes that we do not include, the power law is 
nearly the same but the coefficient is different in the relation between the mean density and 
large separation for different modes. 

We hnd the relation between Az/ and y/p for increasing rotation is very similar to the 


relation found by Suarez et al. (2014) and Garcia Hernandez et al. (2015). In all cases, 
the slope is very similar to that of the classical period-root mean density relation. The 
change with mass for the ZAMS models is not as steep, suggesting that the mean density is 
primarily an indication of evolution stage for models in this range. However, the correlation 
may not be unique when the results of a rapidly rotating model are compared with those 
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Figure 7. Relation between Az/ and y/p for the rotating models considered. Also plotted are the relationships found by previous 
studies for comparison. The purple line shows how this relationship looks if the non-rotating model with 2Mq is evolved and 
the star symbol along this line shows the mean density of a 2Mq ZAMS model rotating at critical rotation. A colour version 
of this figure can be found online. 


of a slightly evolved non-rotating star, as is shown in Figure It goes without saying that 
obtaining a useful estimate of the mass depends crucially on the quality of the surface radius 
information. 


We have also examined how individual frequencies relate to the mean density. Lovekin 


et al. (2009) found that a “mean density” computed using values of the surface radius at 


specihc colatitudes gave a nearly linear relationship between the period and mean density 
when the radius dehned at a colatitude of about 40 degrees was used. No particular signif- 
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icance was attached to this result. We have explored this further and agree with Lovekin 


et ah (2009) for the mode they used, but hud that the colatitude required depends on the 


specihc mode being examined. Specihcally, we hud that either larger n or larger ^ requires 
the colatitude for the surface radius to be used in the “mean density” to increase (i.e., the 
radius must increase). We attach no physical signihcance to this except as an expression that 
models with constant mass and evolutionary state but changing rotation do not constitute 
a homologous family. 


6 FINAL CONSIDERATIONS 


We have examined the variation of both g mode and p mode oscillation frequencies in 
intermediate mass, uniformly rotating, ZAMS stellar models. The behavior of the p mode 
frequencies with rotation is similar to that found for a IOq model by Lovekin & Deupree 


(2008), suggesting that the behavior will be similar for models with masses in between. 


By comparing models of different masses with the same surface shape, we found that the 
period-root mean density relation applies to the p modes for rotating models to the same 
extent and with the same caveats as it does to non-rotating models. 

The g mode frequencies are much less affected by rotation because their frequencies are 
determined in the region of the model between the convective core boundary and the lower 
boundary of the second helium ionization region. The region closer to the convective core is 
probably more important, but the entire region contributes to the frequency determination. 
For sufficiently low frequencies, the frequency ratio of a given mode for two masses is almost 
independent of the mode. For these modes there is a slight trend of increasing frequency 
ratio with decreasing frequency as shown in Figure]^ 

The most complex region of the frequency spectrum is the region of overlap of the p and 
g modes. Because the g mode frequencies increase with rotation and the p mode frequencies 
decrease, the possibility of avoided crossings is high and the frequency ratio for specihc 
modes between two masses may vary signihcantly as a function of rotation rate. Figure 
suggests that modes with higher latitudinal variation show this feature more strongly than 
modes with lower latitudinal variation. 

Another situation in which the scaling with mass will fail is when a given mode at one 
mass has undergone an avoided crossing while it has not at the other mass. If the modes 
rather than the shape of the latitudinal variation can be identihed, the scaling holds. Once 
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the mode has undergone the avoided crossings for both masses, the scaling is restored. 
Despite these problems, the scaling relationships found here should give general guidance 
on how to alter a model to obtain a better approximation to an observed mode frequency 
spectrum. 

There remains the issue of matching modes observed in stars to modes computed for 
specihc models. We have shown that there are certain scaling relations for the frequencies 
of a given mode from one model to another which may be helpful as part of this process. 
However, there are many gaps to £11 before we can routinely associate a given observed mode 
with a specific computed mode. First and foremost is the need to have some reasonable 
estimate of the inclination of the rotation axis to the observer. If the inclination is low, one 


might use the shape of the wings of line profiles as for Vega (Gulliver et ah 1991 Takeda et ah 


2008 Yoon et ah]|2010 ). Stars with very high observed values of V sin i are presumably seen 


at high inclinations. Best is the determination of the inclination by optical interferometry, 
so far restricted to relatively few stars. Even with a reasonably determined inclination, the 
task is still formidable. An examination of the echelle diagram comparing observed mode 


frequencies with computed mode frequencies (e.g.. Figure 5 of Deupree (2011b) for a Oph) 
provides some indication of how far we have yet to go. Clearly, the density of computed modes 
requires a windowing process to help us associate an observed mode with a computed mode. 
For stars with a sizeable number of observed modes, patterns in the observed frequencies 
may exist for comparison with patterns in the computed frequency spectrum, much as for 
the Sun. Stars with only a very limited number of modes with sizeable amplitudes may 
have variations in the colour difference between maximum and minimum light in different 
filters or variations in the line profiles that could help constrain the possibilities. Similarly, 
theoretical determination of which modes are pulsationally unstable or otherwise of such 
low amplitude as to be unobservable at the appropriate inclination would be helpful. 
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